use "$path\datasets\rd_dataset.dta", clear

local depvar "ch_total_child_places_pc"




bandwidth_and_weights, depvar(`depvar') var(margin_1)  bwmethod(CCT) kernel(tri) degree(1) number_of_elec(num_of_obs)




ivreg2 `depvar' elected_women margin_1 inter_1 csu spd greens   if abs(margin_1)<$bw_opt  [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m1
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)
gen esample1=e(sample)






ivreg2 `depvar' elected_women margin_1 inter_1 higher_degree   if abs(margin_1)<$bw_opt  [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m2
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)
gen esample2=e(sample)



ivreg2 `depvar' elected_women margin_1 inter_1  selfemployed student retired    if abs(margin_1)<$bw_opt  [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m3
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)
gen esample3=e(sample)




ivreg2 `depvar' elected_women margin_1 inter_1  alter   if abs(margin_1)<$bw_opt  [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m4
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)
gen esample4=e(sample)

***************************************************************************************************************************


ivreg2 `depvar' elected_women margin_1 inter_1    if abs(margin_1)<$bw_opt & esample1==1  [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m5
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)







ivreg2 `depvar' elected_women margin_1 inter_1    if abs(margin_1)<$bw_opt & esample2==1  [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m6
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)




ivreg2 `depvar' elected_women margin_1 inter_1     if abs(margin_1)<$bw_opt  & esample3==1 [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m7
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)





ivreg2 `depvar' elected_women margin_1 inter_1     if abs(margin_1)<$bw_opt & esample4==1 [pw=weight] , r cluster(gkz ) partial(margin_1 inter_1 )
est store m8
estadd local bw "CCT"
estadd local degree "Linear"
estadd local bw_length  $bw_opt
unique gkz_legis_party if e(sample)
estadd local num_of_elections  `"`r(sum)'"'
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)

************************************************************************************************************************************

esttab  m1 m2 m3 m4 using TableA8/tableA8_PanelA.txt, style(tab) replace order( ) mlabel(,none) ///
cells(b(label(coef.) star fmt(%8.3f) ) se(label((z)) par fmt(%6.3f))) ///
collabels(none) ///
keep (elected_women csu spd greens higher_degree selfemployed student retired   alter  ) ///
stats(bw bw_length degree N   N_clust mean_depvar sd_depvar , layout( @ @ @ @ @  `""@ (@)""' )  fmt( %~#s %9.2f %~# %9.0g %9.0g  %9.2f %9.2f  ) ///
labels("Bandwidth type" "Bandwidth size" "Polynomial"   "N"  "Municipalities" "Mean (SD)"  )) ///
starlevels(* 0.10 ** 0.05 *** 0.01) ///
varlabels( elected_women "Female victory" csu "CSU" spd "SPD" greens "Greens" higher_degree "Higher degree" employed "Employed" selfemployed "Self-employed" student "Student" retired "Retired" housewifehusband "Housewife/-husband" alter "Age" employed_in_any_form "(Self-) employed" ) 



************************************************************************************************************************************

esttab  m5 m6 m7 m8 using TableA8/tableA8_PanelB.txt, style(tab) replace order( ) mlabel(,none) ///
cells(b(label(coef.) star fmt(%8.3f) ) se(label((z)) par fmt(%6.3f))) ///
collabels(none) ///
keep (elected_women   ) ///
stats(bw bw_length degree N   N_clust mean_depvar sd_depvar , layout( @ @ @ @ @  `""@ (@)""' )  fmt( %~#s %9.2f %~# %9.0g %9.0g  %9.2f %9.2f  ) ///
labels("Bandwidth type" "Bandwidth size" "Polynomial"   "N"  "Municipalities" "Mean (SD)"  )) ///
starlevels(* 0.10 ** 0.05 *** 0.01) ///
varlabels( elected_women "Female victory" csu "CSU" spd "SPD" greens "Greens" higher_degree "Higher degree" employed "Employed" selfemployed "Self-employed" student "Student" retired "Retired" housewifehusband "Housewife/-husband" alter "Age" employed_in_any_form "(Self-) employed" ) 
************************************************************************************************************************************
